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A variational theory is proposed for the quasiparticle excitation in high T c cuprates. The theory 
goes beyond the usual Gutzwiller projected mean field state description by including the spin- 
charge recombination effect in the RVB background. The spin-charge recombination effect is found 
to qualitatively alter the behavior of the quasiparticle charge as a function of doping and cause 
considerable anisotropy in quasiparticle weight on the Fermi surface. 



PACS numbers: 

I. INTRODUCTION 

High temperature superconductors are doped Mott in- 
sulators. For such a strongly correlated electron system, 
it is of fundamental importance to know if and how does 
the Landau quasiparticle excitations emerge and what is 
peculiar about their properties in case of their existence. 
In recent decades, experiments, especially the angle re- 
solved photoemission(ARPES) and the scanning tunnel- 
ing microscopic spectrum(STM) measurements have elu- 
cidated in great details the evolution of the quasiparticle 
properties of the high-T c cupartes as a function of doping 
and temperature However, a comprehensive under- 
standing of these experimental observations is still in its 
infancy stage. 

The existence of the quasiparticle excitation is now 
well established in the superconducting state of the 
cuprates [lj . In the overdoped regime, the observations 
basically agrees with our expectations for a conventional 
Landau Fermi liquid, where a large Fermi surface enclos- 
ing a volume consistent with the Luttinger theorem for 
a system with 1 — x electron per unit cell is observed. 
Well defined quasiparticle peaks are found all around the 
Fermi surface too. However, below optimal doping, parts 
of the Fermi surface around the antinodal region become 
less and less clear with decreasing doping. This phenom- 
ena, which is dubbed as nodal-antinodal dichotomy, is 
tangled with the issue of the existence of two gap scales 
and small hole pocket Fermi surface 5-7]. While in the 
half filing limit, in which the cuprates develops antifcr- 
romagnetic long range order, it is generally believed that 
quasiparticle excitation exists in the form of spin polaron 
and form small hole pocket Fermi surface. How is this 
small hole Fermi surface evolved into the large electron 
Fermi surface in the overdoped regime is long unresolved 
problem in the study of the high temperature supercon- 
ductivity. 

Another issue about the quasiparticle excitation in 
high Tc cuprates is their electrodynamic response ker- 
nel, or more loosely, the quasiparticle charge 8]. As a 
result of the Mott physics, the zero temperature super- 



fluid density of the cuprates scales roughly linearly with 
the density of the doped holes x, rather than the result 
of 1 — x when the strong correlation effect is neglected. 
Take it literally, this would imply that in the ground 
state each electron carry a charge of order x rather than 
one. So naively one would expect that the quasiparti- 
cle excitation above the ground state should also carry a 
charge of order x. However, measurements of the temper- 
ature dependence of the superfluid density indicates that 
this is not the case. The almost doping independence of 
the slope of the superfluid density curve as a function of 
temperature indicates that the electrodynamic response 
kernel of the nodal quasiparticle is almost a constant in 
the underdoped regime. 

The total mobile charge density of a doped Mott in- 
sulator is given by density of the doped holes. This re- 
quirement can be trivially satisfied by a small hole Fermi 
surface enclosing an area of a;, on which each quasiparti- 
cle carries a charge of order one, as happens in the anti- 
ferromagnetic ordered state. However, according to the 
Landau Fermi liquid theory, a system with an electron 
density of 1 — x per unit cell and no symmetry breaking 
should always form a large Fermi surface enclosing an 
area of 1 — x. So if we neglect the possible anisotropy of 
the quasiparticle charge on the Fermi surface, one would 
be naturally led to the conclusion that the quasiparti- 
cle on the Fermi surface should each carry a charge of 
order x. We are thus in a situation of dilemma. The 
system should either violate the Luttinger theorem by 
possessing a small Fermi surface of area x in the absence 
of symmetry breaking, or, while having a large Fermi 
surface that respects the Luttinger theorem, have strong 
anisotropy on it. Both of these scenarios points to the 
nontrivial nature of the quasiparticle properties in the 
cuprates. It should be noted that suggestions have been 
made that in a system with topological order, the Lut- 
tinger theorem should be modified to accomodate the 
topological degeneracy [t| [Io[ and a small Fermi surface 
in the absence of symmetry breaking is in principle possi- 
ble. However, up to date there is no evidence in support 
of such a exotic order in the high T c cuprates. 
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The key to the physics of a doped Mott insulator is the 
local constraint of no double occupancy between the elec- 
trons. The slave Boson mean field theory, in which such 
a local constraint is treated in an average manner, pre- 
dicts a lar ge F ermi surface consistent with the Luttinger 
theorem I 111). The quasiparticle weight and quasiparti- 
cle charge is however isotropic on the Fermi surface and 
both show linear scaling with the hole density x. At- 
tempts has been made to go beyond the mean field the- 
ory by including the fluctuation effect but the answer is 
still inconclusive Q. The fluctuation around the mean 
field configuration, which takes the form of gauge fluctu- 
ations, is notoriously hard to be analyzed as there is no 
mass term to control their behavior. 

The origin of the gauge degree of freedom in the slave 
Boson formalism is just the no double occupancy con- 
straint between the electrons. To account for such a lo- 
cal constraint, Gutzwiller projected mean field state of 
the form PG7k Jd— BCS) has been studied extensively 
to address the quasiparticle problem flil - Eoj ]. However, it 
is found that the Gutzwiller projected state inherits to a 
large extent the properties of the mean field theory. For 
example, it also predicts a large Fermi surface on which 
the quasiparticle weight is uniform. Both the quasiparti- 
cle weight and the quasiparticle charge still vanish in the 
half filling limit, albeit with some powers different from 
the mean field theory predictions. 

Despite these problems, both the slave Boson mean 
field theory and the Gutzwiller projection scheme do cap- 
ture one key feature of the cuprates as a doped Mott in- 
sulator: the strong particle-hole asymmetry in the vicin- 
ity of the Fermi energy, which is hard to envisage in 
alternative scenarios where strong correlation effect is 
ignored [ll. 13 . 20j . For this reason, we believe that the 
slave Boson mean field theory and the Gutzwiller pro- 
jected state are good starting point for a more consistent 
theory. 

An important consequence of the gauge fluctuation be- 
yond the mean field description is the so called spin- 
charge recombination effect [8|, [21 1. The gauge fluctua- 



tion, which couples equally to the spinon and the holon 
degree of freedom in the slave Boson formalism, will in- 
duce mutual backflow effect between the two parts. Such 
a backflow effect, which manifests itself as a kinematic 
effect of the local constraint, is beyond the reach of the 
Gutzwiller projection description and encourages the for- 
mation of spinon-holon bound state, whose motion is less 
affected by the gauge fluctuation. 

To account for such a spin-charge recombination effect, 
a RPA theory has been proposed previously [2^ . The the- 
ory introduced a phenomenological attractive interaction 
between the spinon and the holon to induce spinon-holon 
bound state which is interpreted as the quasiparticle exci- 
tation of the system. Although the theory does have the 
potential to explain certain features of the experiments, 
it is not clear to what extent the predictions made by it 



are gauge invariant. The prediction power of the theory 
is also limited by its phenomenological nature and the 
sum rule for the electron spectral function is in general 
violated. Alternatively, a dopon-spinon formalism is in- 
troduced to account for the spin-charge recombination 
effect in which the bare hole is treated as an elementary 
degree of freedom 231] . However, variational calculation 
based on this formalism can only be proceeded in the 
half filling limit with a small and finite number of doped 
holes. 

In this paper, we extend the Gutzwiller projection de- 
scription of the quasiparticle by including the backflow 
effect between the spinon and the holon at the wave 
function level. Our approach has the advantage that its 
predictions are gauge invariant and the electron spec- 
tral function calculated from it satisfy the relevant sum 
rule. We found the quasiparticle charge after the back- 
flow effect correction approaches to a constant value in 
the half filling limit, rather than vanishes as predicted 
by the mean field theory and the Gutzwiller projection 
scheme. At the same time, the quasiparticle weight is 
found to show more and more strong anisotropy on the 
Fermi surface with decreasing doping and by fine tun- 
ing of the Hamiltonian parameters one can indeed repro- 
duce the nodal-antinodal dichotomy phenomena in the 
cuprates. 

This paper is organized as follows. In the next section, 
we review the variational description of the quasiparticle 
excitation in the Gutzwiller projection scheme and moti- 
vate our new variational wave function by reformulating 
the Gutzwiller projected wave function in terms of the 
slave Boson language in which the spin-charge recombi- 
nation effect can be easily included. In section III, we 
introduce the numerical algorithm to do calculation on 
the wave function we have proposed. In section IV, we 
present our numerical results on both the doping depen- 
dence of the quasiparticle charge and the anisotropy of 
the quasiparticle weight on the Fermi surface. In sec- 
tion V, we present some discussion on the results. The 
appendix contains some details of the derivations in the 
text. 



II. VARIATIONAL QUASIPARTICLE WAVE 
FUNCTION 

The model studied in this paper is the t — J model. 

1 



H 



(1) 



in which Cj )(T is the electron operator satisfying the no 
double occupancy constraint J2a ^ia^i,^ — 1> denotes 
the hopping integral between site i and site j. Ylu j\ de- 
notes sum over nearest neighboring sites. In this study, 
ti j will be assumed to be nonzero only between near- 
est neighboring, next nearest neighboring and next next 
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nearest neighboring sites. The corresponding hopping 
integral will be denoted as t, t and t . 

The Gutzwiller projected BCS mean field state of the 
form 



|*) =P A r e P G |BCS), 



(2) 



is generally believed to be a good starting point for a vari- 
ational description of the ground state of the t-J model. 
Here Pg denotes the Gutzwiller projection into the sub- 
space satisfying the constraint ^Z a c\ a Ci^ < 1 and Pjv e 
denotes the projection into the subspace of N e electrons. 
|BCS) is the usual BCS mean field ground state with 
d-wave pairing. 

A variational description of the quasiparticle excita- 
tion on |W) can be constructed in the same spirit as \%f) 
by Gutzwiller projection of the mean field excited state. 
For example, the variational wave function for the quasi- 
particle excitation of hole type has the form, 



|k,a) 



JV.-l 



P G7 ^|BCS), 



in which 7^ a = u^cl a + v^c-^^. This state has been 
studied by many authors in recent years. It is shown that 
|k, a) inherits to a large extent the properties of the mean 
field excitation. In particular, the quasiparticle weight is 
predicted to be isotropic on the underlying Fermi surface. 
At the same time, both the quasiparticle weight and the 
quasiparticle charge response kernel vanish in the half 
filling limit, albeit with some powers different from that 
predicted by the mean field theory. 

To go beyond the Gutzwiller projected wave function 
and to make connections with the effective field theory 
considerations, we reformulate the Gutzwiller projection 
in the slave Boson language. In the slave Boson formu- 
lation, the electron operator is expressed in terms of the 
Fcrmionic spinon operator and the Bosonic holon opera- 
tor as 



(4) 



The no double occupancy constraint now takes the form 
of an equality 



(5) 



In the slave Boson language, the t-J model takes the 
form 



H = -52t id (fl„f jtir btb i + h.c.) 



+ 9 ^""^ ( fi, a fj, a'fj,v' fi , & fi,crfi,<rfj,<j'fj,v')- (6) 

In the mean field treatment, the ground state of the t — J 
model is given by the product of the BCS mean field state 



for the Fermionic spinon and the Bose-Einstein conden- 
sate of the holon. The variational ground state Y&) can 
be shown to be given by the Gutzwiller projection of such 
a product state. 



I*) = P JVh P 8 |/-BCS)®|6-BEC), 



(7) 



in which \f — BCS) is the BCS mean field state of the 
spinon and \b — BEC) is the Bose condensate of the holon. 
For notational convenience, in the following we will ab- 
breviate |/ - BCS) ® \b - BEC) as |SBMF). Here P s 
denotes the projection into the subspace satisfying the 
constraint J2a fi <rfi,<? + b\bi = 1 and Pjv h denotes the 
projection into the subspace with Nh doped holes. The 
BCS state for spinon is of the form 



i/-bcs)= h>' ./>/:,/•;;. 



(3) in which a(i — j) 



(N-N h )/2 



- R i 



(8) 



Here ^ = 
— V tV P *k-(Ri-Rj) _ „v A _ 



^Ei,j A i,j e<k ' (Ri ~ R ^- l lo and A li arc ho PP in g and 
pairing parameters determining the mean field ground 

state and are treated as variational parameters to be opti- 
mized from the variational energy, /j, v is also a variational 
parameter and not to be mistaken as the real chemical 
potential. In this work, ■ will be assumed to have the 
same range as the real hopping integral Uj and A" • will 
be assumed to take the standard d-wave form. 

In the slave Boson language, an electron becomes a 
composite object. To create a hole in the system, one 
should generate a Bogliubov quasiparticle in the BCS 
mean field ground state of the spinon and at the same 
time add a holon to the system. The added holon can 
either enter the Bose condensate of the holon or stay out 
of it. The first choice for the holon leads to the coherent 
quasiparticle peak in the electron spectral function. The 
variational wave function |k, a) for the quasiparticle is 
just given by the Gutzwiller projection of this mean field 
state, i.e., 



|k,a) 



N h + lf s 



7k >Ul SBMF >< ( 9 ) 



here 7k, CT = U *fk,a + w k/-k,ff, 

Up to this point, the slave Boson language seems to 
generate no new result beyond the usual Gutzwiller pro- 
jection scheme. To see the key difference between the 
two schemes, we note that in the usual Gutzwiller pro- 
jection scheme, the commutator between the electron op- 
erator and the Gutzwiller projection operator is nonzero, 
[Pg, Ci.a] 7^ 0, while in the slave Boson language, the 
commutator between the electron operator and the pro- 
jection operator P s is identically zero as a result of the 
gauge invariance of the electron operator. Such a prop- 
erty can be very useful when proving certain sum rules 
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that will be exemplified below. A simple application of 
this property leads readily to the conclusion that the 
electron spectrum in the particle side is totally coher- 
ent in the Gutzwiller projection scheme(as holons are all 
condensed) . 

One more advantage of the slave Boson language is 
that it provides a bridge between the variational study 
and the effective field theory considerations. For exam- 
ple, to describe the spin-charge recombination effect ar- 
gued in the effective field theory context, we can intro- 
duce the following wave function 

|k,a) scr = P Wh+1 P s ^ q7 t_ qi(T 6t|SBMF}, (10) 
q 

in which <p q can be interpreted as the wave function 
for the relative motion between the spinon and holon. 
We note the form of the wave function is quite general. 
For example, when q = 5 q ,o, the Gutzwiller projected 
state |k, a) is recovered, while when tfi q = i>k-q, the bare 
hole state |k, <t) = Ck jCT P./v e PG|BCS) is recovered(see Ap- 
pendix A). In the following, we will take (f> q as variational 
parameters to be determined by the optimization of en- 
ergy. The spin-charge recombination effect then mani- 
fests itself in the short ranged nature of the optimized 

0q- 

In the effective field theory description, the spin-charge 
recombination effect is argued to be caused by the gauge 
fluctuation, which acts to enforce the local constraint be- 
tween the spinon and the holon degree of freedom. In the 
Gutzwiller projection scheme, such a local constraint is 
enforced a posteriori. By so doing, the kinematic effect 
of the constraint, namely the backflow effect between the 
spinon and the holon motion is totally missed. To make 
this point more clearly, we note the kinetic part of the 
Hamiltonian describes a correlated motion of the spinon 
and holon 

Hk = -£ UMtJi^h + h.c), (11) 

and the spinon current is exactly compensated by the 
backflow current of the holon in each hopping steps. The 
main theme of the present work is to elucidate the correc- 
tion induced by the backflow effect on the quasiparticle 
excitations. 

In the slave Boson language, a hole-like quasiparticle 
is composed of a spinon-holon pair. The backflow effect 
cause momentum transfer between these two parts. In 
momentum space, the kinetic energy part of the Hamil- 
tonian is of the form 

Hk = ^ E [ < ( k + q-p)/k+q, CT A^_A + /i.c.] , 

k.p,q,tr 

(12) 

in which t{k) = ± J2i. 3 ■ tij-e jk ' (Ri-R A The scatter- 
ing between the spinon and holon induced by this 



term will in general lead to state state of the form 
S q 0q7k_ q CT &q|SBMF) This is the reason that motivates 
the variational wave function Eq.(10). 

Two things should be noted here. First, in addition 
to causing scattering between the existing spinon and 
holon pair, the backflow effect can also generate extra 
pairs of spinon and holon from the mean field ground 
state. This can be interpreted as a renormalization of 
the ground state, which is not considered in the present 
work. Such a renormalization effect can in principle be 
taken into account in the variational wave function for 
the ground state to arrive at a more consistent descrip- 
tion of the quasiparticle excitation. Second, it can be 
seen that the momentum dependence of the backflow ef- 
fect depends on the detailed form of the hopping integral 
in the Hamiltonian. Thus some of the results presented 
below are not generic, but depends on the Hamiltonian 
parameters. This is especially the case for the nodal- 
antinodal dichotomy phenomena. However, the quasi- 
particle charge around the nodal point (from where the 
contribution to the in-plane transport is the largest) is 
to a large extent not sensitive to the fine tuning of the 
Hamiltonian parameters. 

III. NUMERICAL ALGORITHMS 

The variational wave function proposed above is com- 
posed of N Slater determinants and can be written as 

|k,er) scr = ^0 q |k,q,<r), (13) 
q 

in which |k,q,a) = P^ +1 P s7 ^_ qiff 6 q |SBMF) form a set 
of strongly correlated basis functions. Unlike the mean 
field states 7 ^_ q CT 6 q |SBMF), |k,q, cr) are no longer or- 
thogonal to each other. Furthermore, it can be shown 
that not all |k, q, a) are linearly independent. The proof 
of this point is left to the appendix. 

In terms of this set of strongly correlated basis func- 
tions, the variational energy for the quasiparticle is given 
by 

Ek _ Sq.q' ^q^q' ( k > ^ PjHt-jjk, q' ; g) 

E q , q '^q'( k >cb°1 k >q',0-) 

The minimization of this expression can be casted into a 
generalized eigenvalue problem of the form [24] 

]TH q , q ^ q = A]rO q , q ^ q , (15) 
q q 

in which H q , q > = (k, q, <7|H t _j|k, q', a) and O q , q ' = 
(k, q, <j|k, q', a). The optimized variational energy is 
given by the lowest eigenvalue \ m in of the generalized 
eigenvalue problem. The optimized variational param- 
eters <f) q is given by the eigenvectors corresponding to 
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More generally, the calculation can be interpreted as 
diagonalization of the Hamiltonian in the set of the 
strongly correlated basis functions |k, q, a). Thus if we 
assume approximate completeness of the basis, we can 
even calculate the full electron spectral function A(k, u) 
as well. The expression for A(k, w) in this approximation 
is given by 



A(k,w) 



^O^rS^-iK-Eg)), (16) 



n,q,q' 



in which <f>™ denotes the normalized eigenvector of 
Eq.(15) with eigenvalue A„. denotes the vector corre- 
sponding to a bare hole on the RVB background and is 
given by <jP^ = i>k-qj Eg denotes the variational ground 
state energy. The derivation of Eq.(16) is given in the 
appendix, in which it is also shown that the spectral 
function so calculated satisfies the sum rule of the form 
J dojA(k, oj) = nk- The existence of such a sum rule par- 
tially justifies the approximate completeness of the basis 
functions |k, q, a). 

The most time consuming part of the present calcu- 
lation is the determination of the Hamiltonian matrix 
elements H qq < and the overlap matrix elements O q! q<. In 



Ref [24] , a highly efficient reweighting technique to reach 
this goal is proposed based on the mutual similarity of 
the basis functions. The algorithm has been discussed in 
details in [24j . Here we will only give a brief overview of 
it. 

To do variational Monte Carlo simulation, we expand 
the basis function in a local basis IRA as 



|k,q,0)=$^Vq(fli)|fl< 



(17) 



In our calculation, we have made a particle-hole trans- 
formation on the down spin electron so that the wave 
function ijj q (Ri) takes the form of a product of a Slater 
determinant from the spinon part and a plane wave from 
the holon part. 

To calculate the overlap matrix elements, one can sim- 
ulate the following expression by Monte Carlo method 



E* 1^)1 



MM 



o 



q,q 



E* \M R i 



(18) 



However, there are of order N 2 such terms to be calcu- 
lated and a direct calculation is very time consuming. In 
24], it is shown that a more efficient and statistically 
more stable way to calculate the overlap matrix elements 
is to simulate the following expression 



0„ 



E, W(Ri 



W(Ri) 



Eq Oq,q 



(19) 



in which W(Ri) = Eq IV'qC-^OI 2 - The most important 
advantage of Eq.(19) over Eq.(18) is that the simulation 



of the all N 2 matrix elements O qiq ' can now be done 
in a single run of the Monte Carlo procedure. Another 
advantage is that the statistical error of Eq.(18) caused 
by the nodes in i/j q (Ri) is now reduced. 

To simulate the expression Eq.(19), we choose an ar- 
bitrary basis function |k, qo,c) as a reference state, then 



W(R i ) = \% (R i )\ 2 Yl 



(20) 



and 



W{Ri) 



Eq 



MM 



(21) 



The Hamiltonian matrix elements can be simulated in 
a similar manner. As Eq.(19), we have 



H„ 



! ^(fl.)xH^(«.) 
' W(Ri) 



EqOq,q EiW(Ri) 

in which H^iZi) - Ej(^t-.j\Rj)MRj)> and 

( MM A" 



(22) 



HV(fl,) 



W(Ri) 



(23) 



Thus to simulate the overlap matrix elements and the 
Hamiltonian matrix elements, we only need to calculate 



the ratios 



MM 
4*10 (M 



and ^ q (^) in each Monte Carlo steps. 



This can be easily done with the inverse update trick as 
|k, q, a) and |k, qo, a) differs with each other by at most a 
pair of spinon and holon states. In addition, these two set 
of ratios form two vectors which make their calculation 
highly parallelized. 

The whole computational procedure can then be sum- 
marized as follows. First, we determine the variational 
parameters in the ground state by optimizing the ground 
state energy. Then we calculate the overlap matrix and 
the Hamiltonian matrix with the reweighting technique. 
With these matrixes, we solve the generalized eigenvalue 
problem. The eigenvector corresponding to the lowest 
eigenvalue is the wanted quasiparticle excitation in this 
scheme. From this wave function we can then calculate 
the quasiparticle properties such as its weight, charge, 
and dispersion relation. We can also use the basis as a 
pseudo-complete one to approximate the electron spec- 
tral function to gain an understanding of the incoherent 
part of the spectral function. 



IV. NUMERICAL RESULTS 



Our calculation is done on a v TV x y N square lat- 
tice with periodic-antiperiodic boundary condition. The 
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Hamiltonian parameters are chosen as follows. In our dis- 
cussion of the nodal quasiparticle properties, which are 
insensitive to the fine tuning of the Hamiltonian parame- 
ters, we assume t = and set j = h, y = —0.25. How- 
ever, when discussing the antinodal quasiparticle proper- 
ties, which are sensitive to the value of t , we will present 
results for both t =0 and t = 0.2t (which is more real- 
istic for the cuprates) . The doping concentration studied 
in this work ranges from 2% to 22%. 



A. The quasiparticle peak and the electron spectral 
function 

To illustrate the spin-charge recombination effect on 
the quasiparticle properties, in Fig.l we plot the eigen- 
values of Eq.(15) in ascending order for a system with 
18 x 18 sites and 36 doped holes. The momentum of the 
quasiparticle is chosen along the nodal direction and just 
below the Fermi surface. The eigenvalues split into an 
isolated pole and a continuum with a finite gap between 
them. The existence of the gap implies the formation of 
spinon-holon bound state. The emergence of the bound 
state can also be seen directly from the Fourier transform 
of <p q , which decreases exponentially with the separation 
between the spinon and the holon. 




50 100 150 200 250 300 350 

n 



FIG. 1: The eigenvalues of Eq.(15) in ascending order for a 
system with 18 x 18 sites and 36 doped holes. The momentum 
of the quasiparticle is chosen along the nodal direction and 
slightly below the Fermi surface, t = is assumed in the 
calculation. 

The quasiparticle peak in the electron spectral function 
is just contributed by this spinon-holon bound state. An 
approximate electronic spectral function at the given mo- 
mentum is shown in Fig. 2, showing clearly the emergence 
of the quasiparticle peak out of incoherent background. 
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FIG. 2: The approximate electron spectral function calculated 
from Eq.(16). The parameters used here are the same as those 
used in Fig.l. The 5-function peaks in Eq.(16) are broadened 
with a width of J in the calculation. 

B. Quasiparticle charge 

The quasiparticle charge can be defined through the 
current carried by it. In general, the current carried by a 
quasiparticle can be written in the form j = qv^, where 
ilk = VkEk is the group velocity of the quasiparticle cal- 
culated from its dispersion relation and q is the charge 
of the quasiparticle. As the nodal quasiparticle has the 
largest velocity in the Brillouin zone and thus dominates 
the in-plane electromagnetic response of the system, our 
discussion of the quasiparticle charge will be restricted 
to the nodal quasiparticles. As the velocity of the nodal 
quasiparticle is almost independent of the hole concen- 
tration in the high T c cuprates, we can use the current 
carried by the quasiparticle as a measure of its charge. 

Since the property of the nodal quasiparticle is insen- 
sitive to the fine tuning of the Hamiltonian parameters, 
we will assume t = in the following calculations. The 
electromagnetic current operator of the t-J model in this 
assumption is given by 

jx = it'Yl^l^Ci+x^ - h.c.) 

i,<7 

+ tf^iclaCi+x+y,* - h.C.) 
i.<j 

+ it'^iclad+x-y^-h.c.), (24) 

i.er 

with j y given by a similar expression. In the slave Boson 
mean field theory, the current can be easily found to be 
given by j — xVkek, where ek is the mean field dispersion 
of the quasiparticle. Thus in the slave Boson mean field 
theory, the quasiparticle charge is proportional to the 
hole density x. 



7 



- Spinon-holon recombined 

- Gutzwiller Wave Function 




0.10 0.15 
doping 



FIG. 3: The quasiparticle charge of the nodal quasiparticle as 
a function of the hole concentration. The calculation is done 
on a 14 x 14 lattice, t — is assumed in the calculation. 



The quasiparticle current calculated from our varia- 
tional wave function is shown in Fig. 3 in which the re- 
sult is compared with the predictions of the slave Boson 
mean field theory and the Gutzwiller projected varia- 
tional wave function |k, a). Unlike prediction of the mean 
field theory and the Gutzwiller projection scheme, the 
quasiparticle charge calculated from our wave function 
approaches to a finite value in the half filling limit. A 
non- vanishing quasiparticle charge, which is consistent 
with experimental observation in underdoped cuprates, 
constitutes the main achievement of the present theory. 



weight in the present theory reads 
|o(k,o-|k,o% 



Zy = 



(25) 



As the properties of the off-nodal quasiparticles are 
sensitive to the Hamiltonian parameters, we will present 
results calculated for both the t — 0.2t and t = 
case. The results for t — 0.2t is shown in Fig. 5. Un- 
like the Gutzwiller projected wave function, our varia- 
tional wave function predicts a quasiparticle weight that 
peaks on the Fermi surface. At the same time, the quasi- 
particle weight is anisotropic on the Fermi surface. The 
anisotropy is found to increases with decreasing doping 
and at x = 6% the quasiparticle weight in the nodal re- 
gion is almost three times larger than that in the antin- 
odal region. These predictions of our theory resemble 
closely the experimental observations. 

However, the anisotropy of the quasiparticle weight on 
the Fermi surface is not a generic property of our the- 
ory, but depends on fine tuning of Hamiltonian parame- 
ters. In Fig. 6 we show the quasiparticle weight calculated 
with t =0. The anisotropy of the quasiparticle weight 
on the Fermi surface is found to be much smaller than 
that calculated with t = 0.2t. Thus, in our theory the 
nodal-antinodal dichotomy is not a generic consequence 
of the spin-charge recombination effect, but depends on 
fine tuning of Hamiltonian parameters. The same conclu- 
sion is also reached by the variational calculation based 
on the dopon-spinon formalism (2^. 



V. DISCUSSIONS 



C. The quasiparticle weight and the nodal-antinodal 
dichotomy 

In the slave Boson mean field theory, the quasiparti- 
cle weight is given by Zy = xv\. Thus, Z^ is isotropic 
on the underlying Fermi surface and increases with the 
excitation energy below the Fermi surface. Apart from 
some detailed difference in the doping dependence, the 
mean field predictions on the quasiparticle weight are 
to a large extent inherited by the Gutzwiller projected 
wave function |k, a). However, the increase of the quasi- 
particle weight with excitation energy is obviously at 
odds with the experimental observations and our phys- 
ical intuitions. At the same time, measurements have 
detected large anisotropy in the quasiparticle weight on 
the underlying Fermi surface in the underdoped cuprates. 
These two points constitute the main problems with the 
Gutzwiller projected wave function |k, a). 

To see if the backflow effect can cure these problems of 
|k, a), we have calculated the quasiparticle weight from 
|k, cr) scr as a function of momentum. The quasiparticle 



In this work we have studied the consequence of the 
spin-charge recombination effect on the quasiparticle 
properties of the high-T c cuprates. The spin-charge re- 
combination effect, which can be interpreted as a back- 
flow effect between the spinon and the holon degree of 
freedoms, originates from the no double occupancy con- 
straint of the system. We find such a backflow effect will 
induce spinon-holon bound state and will cause qualita- 
tive changes in the quasiparticle properties. 

In the Gutzwiller projected wave function, the no dou- 
ble occupancy constraint is enforced by hand. However, 
such a posteriorly enforced projection can not account 
for the full effect of the local constraint. Especially, the 
kinematic effect of such a constraint on the motion of the 
spinon and holon, namely the backflow effect, is totally 
beyond the reach of such a description. In the effective 
field theory context, the Gutzwiller projection amounts 
to integration of the temporal component of the gauge 
fluctuation with the spatial component of the gauge fluc- 
tuation totally untouched. Such a unbalanced nature 
of the treatment of the gauge fluctuation explains the 
qualitative similarity between the predictions made by 



8 




FIG. 4: The quasiparticle weight in the first quadrant of the Brillouin zone for t = 0.2t calculated from the spin-charge 
recombined wave function((a) and (c)) and the Gutzwiller projected wave function((b) and (d)) at x — 6%((a) and (b)) and 
x — 14% ((c) and (d)). The calculation is done on a 10 x 10 lattice. 



the mean field theory and the Gutzwiller projected wave 
function. 

The most remarkable consequence of the backflow ef- 
fect on the quasiparticle properties is the modification 
of its electrodynamic response kernel. In the mean field 
theory, in which the motion of the spinon and holon is in- 
dependent of each other, a quasiparticle carries a charge 
of x through the holon condensate. The vanishing of 
the quasiparticle charge in the half filling limit predicted 
by the mean field theory is inherited by the Gutzwiller 
projected wave function. After the backflow effect cor- 
rection, the spinon and the holon form bound state. The 
holon dragged by the spinon in such a composite object 
contributes a nonvanishing charge to the quasiparticle 
even in the half filling limit. It should be emphasized 
that this is generic consequence of the backflow effect 
and does not depends on fine tuning of Hamiltonian pa- 
rameters. 

The backflow effect can also provide a potential mecha- 
nism for the experimental observation of nodal-antinodal 
dichotomy and at the same time result in a quasiparticle 
weight that decreases with increasing excitation energy, 
as would be expected from general physical arguments. 



However, before making serious comparisons with exper- 
iments, it should be kept in mind that unlike the nodal 
quasiparticles, the off-nodal quasiparticles are in general 
sensitive to the fine tuning of Hamiltonian parameters. 
In particular, we find the anisotropy of the quasiparti- 
cle weight on the Fermi surface depends crucially on the 
value of t . We thus can not exclude the possibility that 
some other more generic mechanism is responsible for the 
observed nodal-antinodal dichotomy. 

Finally, we note that the backflow effect will also cause 
renormalization of the ground state. In our calculation, 
we have assumed implicitly that such a renormalization 
is not strong enough to induce Fermi surface reconstruc- 
tion, in which case our calculation would be totally in- 
valid. In the absence of the Fermi surface reconstruction, 
such a renormalization effect on the ground state can in 
principle be taken into account in our theory to arrive at 
a more consistent description of the quasiparticle proper- 
ties. When the backflow effect is strong enough to cause 
Fermi surface reconstruction, there arises the interesting 
and exotic possibility of forming a small Fermi pocket 
without any symmetry breaking. An important problem 
then is whether the small Fermi pocket has a quantized 



1.0 0.0 



1.0 0.0 



FIG. 5: The quasiparticle weight in the first quadrant of the Brillouin zone for t — calculated from the spin-charge recombined 
wave function((a) and (c)) and the Gutzwiller projected wave function((b) and (d)) at x — 6%((a) and (b)) and x = 14%((c) 
and (d)). The calculation is done on a 10 x 10 lattice. 



volume. We leave these issues to future study. 

Tao Li is supported by NSFC Grant No. 10774187 
and National Basic Research Program of China No. 
2007CB925001, Fan Yang is grateful for the NSFC Grant 
No.10704008. 

APPDENIX A 

In this appendix, we show the rank of the basis func- 
tions |k, q, a) is lower than that of the mean field states 
7j[_ q |SBMF) by one. This can be shown by proving 
the following identity 

5> k _ q |k,q 7( r) =0, (26) 
q 

in which = |(1 — |r-). 

Using the Bogliubov transformation fl a = Uk7 kcr + 

w k 7_ kjff , we have Uk-q7k-q,a = /k-q,<7 ~ "k-q7-k+ q ,ff- 
As 7-k+ q ,ff|/ — BCS) = 0, we have 

^u k _ q |k,q, a) 
q 



= ^P^+tPs/L^&qlSBMF) 

q 

= ^P^+A/^bJlSBMF) 

i 

= 0, (27) 

in which the no double occupancy constraint has been 
used in the final step. 

APPDENIX B 

In this appendix, we derive the expression for the elec- 
tron spectral function in the approximation that |k, q, a) 
form a complete set for the description of the quasipar- 
ticle excitation and show that the spectral function so 
calculated does satisfy the sum rule. 

A bare hole created in the ground state is given by 
|k, <r)o = Ck l(T P.iv E PG|BCS). Written in terms of the slave 
particles, it reads 

|k, a) = P^+iP s /k-q, ff &q|SBMF) = J2 ^ l k > q> CT X 28 ) 
q q 
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in which 0° = Uk- q . In the derivation we have used the 
fact that P s and fi t ab\ commute with each other. 

Assuming that |k, q, a) form a complete set, the elec- 
tronic spectral function can then be calculated as 

A(k,u) = lo(k, o-|k, a) n \ 2 S{io - (E k ,„ - E fl )), (29) 



in which |k, <j) n denotes the n— th eigenvector of Eq.(15) 
and E k .„ is corresponding eigenvalue. E ff is the varia- 
tional ground state energy. We thus have 

A(k,u) = ]T K(W£| a *(w - (E k ,„ - E s )), (30) 

n,q,q' 

in which </>™ is the n— th eigenvector of Eq.(15) and sat- 
isfies the following orthonormal condition 



q>q' 



(31) 



The electronic spectral function so obtained satisfy the 
sum rule J dojA(k, uj) — (*|t^ <r <Jk lCT |'4 r ). In fact, from 
Eq.(5) we have 



/ duA%u) = lo(k, 

ajk, cr) 



=0 (k, cr|k,cr) 



Ck)ff |*) 



(32) 



where we have used the fact that |k, a) n forms a complete 
set in the space spanned by |k, q, a). 
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